The high-level error bound for shifted 
surface spline interpolation 

Lin-Tian Luh* 

Department of Mathematics, Providence University, 
Shalu, Taichung, Taiwan 
email: ltluh@pu.edu.tw 
phone: (04)26328001 ext. 15126 
fax: (04)26324653 

February 2, 2008 



Radial function interpolation of scattered data is a frequently used method 
for multivariate data fitting. One of the most frequently used radial functions 
is called shifted surface spline, introduced by Dyn, Levin and Rippa in |I] 
for R 2 . Then it's extended to R n for n > 1. Many articles have studied its 
properties, as can be seen in [21 El El HSl HE1 UHl [IHl EI3 EU - When dealing 
with this function, the most commonly used error bounds are the one raised 
by Wu and Schaback in ^7j, and the one raised by Madych and Nelson in 
|13j . Both are 0(d l ) as d — > 0, where I is a positive integer and d is the fill- 
distance. In this paper we present an improved error bound which is 0(u 1 ' d ) 
as d — ► 0, where < uj < 1 is a constant which can be accurately calculated. 
Keywords:radial basis function, shifted surface spline, error bound. 
AMS subject classification: 41A05, 41A15, 41A25, 41A30, 41A63, 65D10. 

1 Introduction 

Let /ibea continuous function on R n which is conditionally positive definite 
of order m. Given data (xj, fj), j = 1, . . . , N, where X = {xi, . . . , xn} is a 
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subset of points in R n and the /-s are real or complex numbers, the so-called 
h spline interpolant of these data is the function s defined by 

N 

s(x) = p(x) +*^2cjh(x — Xj), (1) 
i=i 

where p(x) is a polynomial in P m -\ and c'^s are chosen so that 

JV 

£c,g( % ) = o (2) 

3=1 

for all polynomials q in P m -i and 

N 

p{xi) + Cjh(xi - xj) = fi, i = 1, . . . , N. (3) 
i=i 

Here P m -i denotes the class of those polynomials of R n of degree < m — 1. 

It is well known that the system of equations © and (j3J) has a unique 
solution when X is a determining set for P m -i and h is strictly conditionally 
positive definite. For more details please see [12]. Thus, in this case, the 
interpolant s(x) is well defined. 

We remind the reader that X is said to be a determining set for P m -\ if 
p is in P m -i and p vanishes on X imply that p is identically zero. 

In this paper h is defined by formula 

h(x) := (-l) m (|x| 2 + c 2 )$log(\x\ 2 + c 2 )^, A G Z + , m = 1 + c> 0, 
x e R n , A, n ewen, (4) 

where |x| is the Euclidean norm of x, and A, c are constants. In fact, the def- 
inition of shifted surface spline covers odd dimensions. For odd dimensions, 
it's of the form 

h(x) := (-lY x -^(\x\ 2 + c 2 ) x ~%, nodd, A G Z + = {1,2,3,...} 
n 

ana A > — . (5) 

However, this is just multiquadric whose exponential error estimates have 
already been constructed by Madych and Nelson in (Hj, and the calculation 
of the constant uj in 0(u>d) can be found in JT]. Hence we will not discuss it. 
Instead, we will focus on even dimensions. 
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1.1 A Bound for Multivariate Polynomials 

A key ingredient in the development of our estimates is the following lemma 
which gives a bound on the size of a polynomial on a cube in R n in terms 
of its values on a discrete subset which is scattered in a sufficiently uniform 
manner. We cite it directly from ^3] and omit its proof. 

Lemma 1.1 For n = 1,2,..., define 7„ by the formulae 71 = 2 and, if 
n > 1, j n = 2n(l + 7 n _i). Let Q be a cube in R n that is subdivided into q n 
identical subcubes. Let Y be a set of q n points obtained by selecting a point 
from each of those subcubes. If q > 7„(fc + 1), then for all p in Pf. 

sup x£Q \p{x)\ < e 2n ^ {k+1) sup yeY \p(y)\. 

1.2 A Variational Framework for Interpolation 

The precise statement of our estimates concerning h splines requires a certain 
amount of technical notation and terminology which is identical to that used 
in [T^j. For the convenience of the reader we recall several basic notions. 

The space of complex-valued functions on R n that are compactly sup- 
ported and infinitely different iable is denoted by T>. The Fourier transform 
of a function in T> is 

0(f) = J e- l<x 'Z><f){x)dx. 
A continuous function h is conditionally positive definite of order m if 

J h(x)(f)(x) * (f)(x)dx > 

holds whenever <fi = p(D)ip with tp in T> and p(D) a linear homogeneous 
constant coefficient differential operator of order m. Here = 4>(—x) and * 
denotes the convolution product 

01 * <j>2(t) = J (j>l(x)(f>2(t — X)dx. 

As pointed out in ^3], this definition of conditional positive definiteness is 
equivalent to that of [TJ] which is generally used. 
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If h is a continuous conditionally positive definite function of order m, 
the Fourier transform of h uniquely determines a positive Borel measure /i 
on i? n \{0} and constants a 7 , |7| = 2m as follows: For all ip G T> 

fh{ X )^x)dx = /{^(e)-x(e) E ^(0)^)^(0 

| 7 |<2m 

+ E ^(0)^' ( 6 ) 

where for every choice of complex numbers c a , |ot| = to, 

E E a a+pCefip > 0. 

|o|=m |/3|=m 

Here \ is a function in X> such that 1 — x(^) has a zero of order 2m + 1 at 
£ = 0; both of the integrals j*o<|£|<i |£| 2m c£/i(£), /^|>i^A i (0 are finite. The 
choice of x affects the value of the coefficients a 7 for I7I < 2m. 

Our variational framework for interpolation is supplied by a space we 
denote by Ch, m - If 

"Cm = {<t> £V : J x a (f)(x)dx = for all \a\ < m}, 

then Ch m is the class of those continuous functions / which satisfy 



f(x)(j)(x)dx 



<c(/) h(x-y)<l>(x)<t>{y)dxdy\ (7) 



for some constant c(f) and all in V m . If / G C/i, m , let denote the 

smallest constant c(f) for which (7) is true. Recall that ||/|| is a semi-norm 
and Ch, m is a semi-Hilbert space; in the case m = it is a norm and a Hilbert 
space respectively. 

2 Main Results 

We first recall that the function h defined in (4) is conditionally positive 
definite of order m = 1 + |. This can be found in and many relevant 
papers. Its Fourier transform jB] is 

h(6) =l{\,n)\6\- x - n }Cn±x{c\6\) (8) 
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where l(X, n) > is a constant depending on A and n, and lC u (t) = t u K, u (t), 
K, v (t) being the modified Bessel function of the second kindPQ. Then we have 
the following lemma. 

Lemma 2.1 Let h be as in (4) and m be its order of conditional positive 
definiteness. There exists a positive constant p such that 

I \i\ k d^)<l{\n)- x [^--n-a n -c x - k -^-p k -k\ (9) 

JR n V 2 

for all integer k > 2m + 2 where p is defined in (6), a n denotes the volume 
of the unit ball in R n , c is as in (4), and A is a positive constant. 

Proof. We first transform the integral of the left-hand side of the inequality 
into a simpler form. 

R n 

\Z\H(\,n)jC^(c\m\- X - n dZ 

R n 2 

JR n 2 

n+A / 7T / | j_ | ]L n + A 



i(A,n)cT" t /- / |{|"-— • -^=— d{ 



2Jr 

£(A,n)c 2 ■ \ — ■ n ■ a n / r 2 • — = ctr 

v ' y V 2 Jo x/cr • e cr 



n + A /7f 1 f f+ 2 

l{\n)c * J- ■ n- a n - —= — dr 

V 2 y'c Jo e cr 

t, | ra — A — 3 

Z A, n)c 2 \ --n-a n - —=■ - ■ „_ A-1 / dr 

V 2 ^/c c fc + — 2 — Jo e r 

n W \-k f°° rk> , , ,/ , ra-A-3 

l{ A, n)\ — ■ n ■ a n ■ c / — dr where k =k + 
V 2 Jo e r 
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Note that £;>2m + 2 = 4 + A implies k' > ^±f±^ > 0. 

Now we divide the proof into three cases. Let k" = \k'] which is the 
smallest integer greater than or equal to k! . 

Casel. Assume k" > k. Let k" = k + s. Then 

POO POO y h 

\ —dr< dr = k"\ = (k + s)(k + s — V) ■ • ■ (k + l)k\ 

Jo e r Jo e r 



and 

/ dr< dr = {k" + l)\ = {k + s + l){k + s)---{k + 2){k + l)k\. 

Jo e r Jo e r 

Note that 

(k + s + l)(k + s) ■ ■ ■ (k + 2) _ k + s + 1 
(k + s)(k + s- l)...(fc + 1) ~ k + 1 

The condition k > 2m + 2 implies that 

k+s+1 2m + 3 + s , s 

< = H . 

k + 1 ~ 2m + 3 2m + 3 

Let P = 1 + 2^3- Then 



PCX) +1 

/ ^~dr< A ■ p k+1 ■ (k + 1)\ 
Jo e 



if / °° V rfr < Ao-p fc -A;!. The smallest fc" is = 2m + 2 + s when k = 2m + 2. 
Now, 



-.k r . 



PCX) y^Q 

/ — dr = fcg! = (2m + 2 + s)(2m + 1 + s) • • • (2m + 3) (2m + 2)! 

j o e r 

(2m + 2 + g )(2m + l + s)---(2m + 3) 2m+2 

= A -p 2m+2 -(2m + 2)! 

(2m + 2 + s)(2m + 1 + s) ■ ■ ■ (2m + 3) 



where A 



p2m+2 



It follows that Jo°° r —dr < A ■ p ■ k\ for all k > 2m + 2. 

Case2. Assume k" < k. Let k" — k — s where s > 0. Then 

h' h" 1 

/•OO ip^ /*00 rp h I 

/ —dr< dr = k"\ = (k-s)\ = — r Tl r-k\ 

Jo e r - Jo e r v ; k (k - 1 )■■■ (k-s + l) 



and 

r k'+l 

PCX) 

/ — dr < / 

Jo e r Jo 



-dr 



& 
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Note that 



( (k + l)k ■ ■ ■ (k - s + 2) / - 1) • • • (k - s + 1) 
k(k - 1) • • • (k - s + 1) 
(k + l)k---(k-s + 2) 
(k-a + l) 
k + 1 

< 1. 



Let p = 1. Then 

„fc"+i 



if / °° r —dr< A ■ p k ■ k\. The smallest k is k = 2m + 2. Hence the smallest 
k" is k% = k - s = 2m + 2 - s. Now, 



h" 

POO y^Q 

/ — dr = kZ\ = (2m + 2-s)\ = (ko-s)\ 
Jo e r 



fc (A;o - 1) • • • (fo - s + 1) 
= A • p ko ■ k \ where A = 



(2m + 2) (2m + 1) • • • (2m - s + 3) ' 



Jo " 

Case3. Assume k" = k. Then 



It follows that / °° < A • p k ■ k\ for all k > 2m + 2. 



/ — < / dr = k\ and / dr < (k + 

Jo e r ~ Jo e r Jo e r ~ 

Let p = 1. Then / °° ^dr < A • p k ■ k\ for all k where A = 1. 
The lemma is now an immediate result of the three cases. 



Remark: For the convenience of the reader we should express the constants 
A and p in a clear form. It's easily shown that 

(a) k" > k if and only if n — A > 3, 

(b) k" < k if and only if n — A < 1, and 

(c) k" = k if and only if 1 < n — A < 3, 

where k" and are as in the proof of the lemma. We thus have the following 
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situations. 

, 3 T,p+ « = [So 



a)n - A > 3. Let s = f^^f^] . Then 



s (2m + 2 + s)(2m + l + s)---(2m + 3) 



p = 1 H and A — — — 

H 2m + 3 p 2m + 2 

(b)n - A < 1. Let a = - [^^1 . Then 



p = 1 and A c 



2 

1 



(2m + 2) (2m + 1) • • • (2m - s + 3) ' 
(c)l < n - A < 3. We have 

p = 1 and A = 1. 

Before introducing our main theorem, we need the following two lemmas, 
first of which is cited directly from|14j. 

Lemma 2.2 Let Q, Y, and 7„ be as in Lemmal.l. Then, given a point x 
in Q, there is a measure a supported on Y such that 



p(y)do-(y) =p{x) 



for all p in P k , and 

d\a\(y) < e 2n ^ {k+1) . 



Lemma 2.3 For any positive integer k, 



(2*)! 

<2 k . 



Proof. This inequality holds for k = 1 obviously. We proceed by induction. 



[2(fc + l)]! J(2k + 2)\ J(2k)\ J(2k + 2)(2k + l) 



(Jfc + l)! fc!(Jfe + l) k\ k + 1 

< \/W \/^±5 <2fc (2fc + 2) _ . fc+1 

~ k\ k+1 k + 1 ' H 

Because of the local nature of the result, we first restrict our attention to 
the case where x lies in a cube. 
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Theorem 2.4 Suppose h is defined as in (4)- Let p be its corresponding 
measure as in (6). Then, given a positive number b , there are positive con- 
stants 5q and u), < u < 1, which depend on bo for which the following is 
true: 

If f G Ch, m and s is the h spline that interpolates f on a subset X of R n , 
then 

\f(x) - s(x)\ < sjl(\, n) ■ C^) l,A ■ v^ 7 ^ ■^■^A .J-\\f\\ h (10) 

holds for all x in a cube E provided that (a)E has side b and b > b , (b)0 < 
5 < 5 and (c)every subcube of E of side 5 contains a point of X. Here, 
l(X,n) is defined in (8), a n denotes the volume of the unit ball in R n , and 
c, A are as in (9). 

The number <5 and uj can be expressed specifically as 

1 /2\scfe 

where 

C = max\2p'^e 2n ^, p' = 9 -. 

I 3o J c 

The number p can be found in the remark following Lemma2.1, 7 n is defined 
in Lemmal.l, and m = 1 + | is defined in (4). 

Proof. First, let p, j n , and bo be the constants appearing in Lemma2.1, 
Lemmal.l, and Theorem2.4, repectively. Let 

B = 2p / v / ne 2n7 " and C = max(B, -M 

I 3d J 

where p' = -. Let 

R 1 

3C 7n (m + l)' 

where m is the order of c.p.d. of h. 

Now, let x be any point of the cube E and recall that Theorem4.2 of JH] 
implies that 

|/(aO-a(s)|<cfc||/|| & / \y-x\ k d\a\(y) (11) 
whenever k > m, where a is any measure supported on X such that 

/ p(y)da(y)=p(x) (12) 
9 



for all polynomials p in Pk-i- Here 




whenever k > m. By (9), for all 2k > 2m + 2, 

< y/l(X, n) • . ^ .^■c-'-^A'o- {2p) k (13) 

due to Lemma2.3. 

To obtain the desired bound on \f(x) — s(x)\, it suffices to find a suitable 
bound for 

I = c k f \y - x\ k d\a\(y). 

This is done by choosing the measure a appropriately. We proceed as follows: 
Let 5 be a parameter as in the statement of the theorem. Since 5 < 5q 
and < ?>C^ n 5 < we may choose an integer k > m + 1 so that 

1 < SC-f n k5 < 2. 

Note that 7„£;<5 < b for such a k. Let Q be any cube which contains x, 
has side 7„A;<5, and is contained in E. Subdivide Q into {^f n k) n congruent 
subcubes of side 5. Since each of these subcubes must contain a point of X, 
select a point of X from each such subcube and call the resulting discrete 
set Y . By virtue of Lemma2.2 we may conclude that there is a measure a 
supported on Y which satisfies (12) and enjoys the estimate 

/ d\a\(y) <e 2n ^ k . (14) 

We use this measure in (11) to obtain an estimate on I. 

Using (13), (14), and the fact that support of a is contained in Q whose 
diameter is ^/n^ n k8, we may write 

/ < 0(^-(|) 1/4 -v^-c'-c- fc . ^■(2p) k .(V^lnk5) k e 2 ^ k 
< {C ln k5) k {^l{\n) • (|) x / 4 • ^T n - c * • yfe). (15) 
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Since 

(15) implies that 
I < 



2 1 

C-y n k5 < - and k > — — -, 

3 3C<y n d 



Hence we may conclude that 



\f(x)-s(x)\<^n) 



C2 



where 



a; 



3C 7 n 



This completes the proof. (j 
What's noteworthy is that in Theorem2.4 the parameter 5 is not the 

generally used fill-distance. For easy use we should transform the theorem 

into a statement described by the fill-distance. 
Let 

d(Q,X) = sup ye ^inf xeX \y - x\ 

be the fill-distance. Observe that every cube of side 8 contains a ball of radius 
|. Thus the subcube condition in Theorem2.4 is satisfied when 5 = 2d(E, X). 
More generally, we can easily conclude the following: 



Corollary 2.5 Suppose h is defined as in (4)- Let fi be its corresponding 
measure as in (6). Then, given a positive number b , there are positive con- 
stants do and u', < uj' < 1, which depend on b for which the following is 
true: If f e C^ m and s is the h spline that interpolates f on a subset X of 
R n , then 



\f(x)-s(x)\<^l{X^-(^-^ 



na r 



A 

c 2 



(16) 



holds for all x in a cube E C Q, where Q is a set which can be expressed 
as the union of rotations and translations of a fixed cube of side bo, provided 
that (a)E has side b > bo, (b)0 < d < d and (c)every subcube of E of side 
2d contains a point of X. Here, a n denotes the volume of the unit ball in R n 
and c, A are as in (9). Moreover d = y and uj' = v /cJ where 5q and uj are 
as in Theorem2.4- 
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Proof. Let d = f and 5 = 2d. Then < d < d iff < 5 < S . Our 

ii i i 

corollary follows immediately by noting that w« = ujm = y/u d = (uj')d. (j 

Remark: The space Ch, m probably is unfamiliar to most people. It's in- 
troduced by Madych and Nelson in [T2] and [13]. Later Luh made char- 
acterizations for it in [Zj and [Hj. Some people think that it's defined by 
Gelfand and Shilov's generalized Fourier transform, and is therefore difficult 
to deal with. This is not true. In fact, it can be characterized by Schwartz's 
generalized Fourier transform. The situation is not so bad. Moreover, some 
people think that Ch, m is the closure of Wu and Schaback's function space 
which is defined in ^Tj\. This is also not true. The two spaces have very little 
connection. Luh has also made a clarification for this problem. For further 
details, please see [U] and [TUj . 
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